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Andreev reflection is investigated in layered anisotropic 
normal metal / superconductor (N/S) systems in the case of 
an energy gap (A) in S not negligible with respect to the Fermi 
energy (E_f), as it probably occurs with high critical temper- 
ature superconductors (HFS). We find that in these limits 
retro-reflectivity, which is a fundamental feature of Andreev 
reflection, is broken modifying sensitively transport across 
S/N interfaces. We discuss the consequences for supercur- 
rents in HFS Josephson junctions and for the midgap states 
in S-N contacts. 
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Andreev reflection (AR) is a scattering process occur- 
ring at superconductor/ normal metal (S/N) interfaces 
that converts an electron incident on a superconductor 
into a hole, while a Cooper pair is added to the super- 
conducting condensate.cl Because of conservation of mo- 
mentum, the hole is reflected back in the direction of the 
incoming electron and all components of the velocity are 
substantially inverted if the exchange momentum in the 
scattering process is much less than the Fermi momen- 
tum. Retro-reflection occurs whenever the Fermi energy 
(Ef) is much larger than the gap value (A) (Andreev ap- 
proximation). Such an approximation neglects that the 
retro-reflected hole has in reality a different momentum 
<5k in the direction perpendicular to the_SYN interface, 
which is proportional to the ratio A/i^Era This means 
that retro-reflectivity is broken in some conditions and 
this will be one of the main issues of this paper. 

While predictions based on the Andreev approxima- 
tion provide accurate explanations in systems employing 
low critical temperature superconductors, in high-critical 
temperature superconductor (HTS) structures the sit- 
uation is more questionable. If we consider that the 
gap value could be in some directions of the order of 
20 meV (one order of magnitude larger than A of tra- 
ditional superconductors) and that the Fermi energy is 
roughly one order of magnitude less than of tradi- 
tional superconductors p it is interesting to go beyond 
the Andreev approximation for systems employing HTS. 
Concepts based on AR have been widely used to inter- 
pret propertieS|-pf HTS grain boundary (GB) Josephson 
junctions (JJ).a Some new interesting arguments have 
been developed by taking into account unconventional 



order parameter symmetry. Examples are given by the 
presence of zero bound states in the density of states 
of YBa2Cu307 in the (110) crystaUographic directiora-El 
and more generally by phenomena associated with bro- 
ken time reversal symmetry (BTRS).Q 

In this paper we demonstrate that the effects neglected 
in the Andreev approximation may determine an extreme 
depression of Andreev reflection processes in some di- 
rections at S/N interfaces and an enhancement of or- 
dinary scattering. The implications of these effects on 
bound states at interfaces employing superconductors 
with d-wave order parameter symmetry are also consid- 
ered. These phenomena can reveal several important fea- 
tures in charge transport in HTS JJ and enlighten some 
aspects of the phenomenology of the junctions within 
the framework of fundamental issues of HTS, such as 
anisotropy. We stress that the effect we consider, being 
intrinsically related to Andreev reflection, provides a mi- 
croscopic explanation of an intrinsic enhanced scattering 
at the S/N interface. 

Before taking into account an order parameter with a 
d-wave symmetrya, we consider a layered normal metal 
facing an isotropic superconductor with a high value of 
the order parameter typical for HTS. This is illustrated 
in the junction cross section scheme of Fig. la and in 3- 
dimensional view in Fig. lb. An electron moving along 
the planes tilted at an angle 6\ with respect to the junc- 
tion interface is reflected as a hole at an angle O2 ■ Locally 
Andreev reflection tends to move quasi-particles out of 
plane and to favor in some way transport along the c- 
axis. This counterbalances the fact that quasi-particle 
transport in HTS is much favored along the a-b planes. 
We will present calculations and phenomenological pre- 
dictions for the layered structures reported in Fig.lc and 
Fig. Id. These can be considered representative, of (100) 
and non-ideal (001) tilt GB JJs respectively.ElH 

General formalism: Andreev reflection proba- 
bility for large values of the A/E F ratio. In order 
to describe charge transport through a normal metal - su- 
perconductor (N/S) interface, we have, used the Blonder- 
Tinkham-Klapwijk (BTK) approacbliJ introducing some 
significant modifications. In solving the Bogoliubov - de 
Gennes (BdG) equations for the wave functions ip n ,s m 
the N, S regions, we consider the terms of the order of 
A/Ep in the expressions of the wave vectors along the 
incoming (electronic) and reflected (hole) trajectories in 
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FIG. 1. a) Cross section of a S/N interface with N lay- 
ered normal metal and S isotropic superconductor with a high 
value of the order parameter. An electron moving along the 
planes tilted of an angle 8i with respect to the junction in- 
terface is reflected as a hole at an angle 82 ; b) 3-dim. view of 
(a); c) cross section of a S/N interface like in (a) being both 
electrodes layered, d) Scheme of non-ideal (001) GB JJ along 
with the supercurrent transport mechanism. 



N, k e ^ — kp n y/l ± E/Ep, and along the transmitted 
trajectories without and with branch crossing in S, k± = 



^(0). Here 



k Fs ^ 1 ± Y -E 2 - |A± /Ep (see Fig.la). The expression 

for k± is extended to the general anisotropic case, when 
the magnitude | A ± | and the phase ip± of the gap function 
in S are specified on the trajectories k±. We also gener- 
alize the BTK matching conditions for ip niS at the N/S 

interface: ip n = ip s , ^(ip a - VO = ( H 

Hi, 2 = J (U(x)—E ei h)dx, U(x) is the interface barrier po- 
tential, Eej-t are the kinetic energies of the electrons and 
the holes. In the WKB approximation barrier transmis- 
sion coefficients D e ^ cx exp(— 2 J ^U(x) — E e ,h/frdx) 
and Dh < D e since the— hole has less kinetic energy 
to overcome the barrier S3 In analogy with the BTK 
model, the dimensionless barrier strengths for electrons 
and holes can be introduced Z lj2 = — -D e ,/i)/-D e ji 
and generally Z\ < Z 2 . Shis effect may only hold for a re- 
alistic extended barrierEj rather than for a delta-barrier. 

With these modifications, we solve the BdG equations 
and find the probability of the Andreev and normal re- 
flection process for an electron incoming from the N side 
respectively: 



A(E) = 



B(E) 
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%(l-iy/\A* ± \-EP/E). 

This approach formally presents some analogies with 
the formulation of the problem of spin- polarized tuuacL 
ing in ferromagnet / superconductor (F/S) junctions£3'tj 
A formal analogy can be established for example between 
the Fermi momenta in the spin subbands and the Fermi 
momenta in the planes (fen = k n+ ) and across the planes 
(k± = k n -) respectively. By increasing the mismatch 
between them, in both cases the contribution to the cur- 
rent due to Andreev reflection is reduced. The effects 
considered in the present paper arise from the loss of 
retroreflectivity due to large gap values, while in the F/S 
interfaces they are due to an exchange field in F. 

Josephson current: comparison with HTS JJ. 
The modification of the probability of the Andreev re- 
flection process has direct consequences on the calcula- 
tion of the Josephson current carried by Andreev bound 
states (see Fig. Id). Andreev bound states are localized 
in the barrier region and are formed by an electron and a 
hole moving in opposite directions. As a consequence the 
momenta mismatch between an electron and a hole leads 
to a depression of the Josephson current. As a generic 
case, we consider tunnel SIS junction, Z\ = Z% = Z 3> 1 
and both electrodes as s-wave superconductors. 

Let an electron have an angle a e relative to the planes. 
As discussed in the introduction, the Andreev reflected 
hole moves at an angle ah different from a e . The an- 
gles a e ^h are related by the conservation of the momenta 
parallel to the interface sin(a — a e )k e = sin(a — Qth)kh, 
where a — (tt/2 — 8±) is the angle between the planes 
and the interface normal and the electron (hole) mo- 
menta, k e ,h are given by (1 ± E B /E F )k~l = cos 2 (a e ^ h - 



)fc|| + sin 2 (a e ,h — cv)fcj 



Here Eb is the Andreev 
bound state energy, and the anisotropic Fermi surface is 
approximated by an ellipsoid with axes (k± , feii). The 
results do not depend qualitatively on this choice.|— . 

According to the formalism of Furusaki et al&3 the 
Josephson current per conductance channel is expressed 
via the amplitude of Andreev reflection a(ip,u) n ) in the 
barrier region (N) 



2Ti 



E 



a{(f,u> n ) _ a{-ip,w n ) 
k e kh 



(3) 



Here uj n = nT(2n+ 1) and <p is the phase difference. The 
amplitude a(ip,u> n ) is found from the solution of BdG 
equations and describes the multiple scattering in the 
barrier region. In a tunnel junction Eb = A. We neglect 
here a weak energy dependence of A due to non-constant 
density of states in the relevant energy range. 

The angle dependence of a(ip, ui n ) in Eq.(||) is mainly 
controlled by the scattering amplitude t e h in the electron- 
hole channel, a cx t e h(cx,a e ,ah)- The normal state 
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FIG. 3. The conductance for different values of the A/Ef 
ratio in a tunnel NIS junction. Inset: resonant state origi- 
nating at the interface of a d-wave S facing an insulator or a 
normal metal N' along the (110) orientation. 

tion interfacesflEiEIll deserves further investigations. 
Interface roughness produces a smooth dependence of Ic 
on the angle a, in contrast to the sharp dependences 
shown above. 

S (anisotropic superconductor) / N interface: 
conductance and the problem of zero bound state. 

We consider another aspect typical of HTS junctions re- 
lated to Andreev reflection by introducing a d-wave or- 
der parameter for the S and investigating the origin of 
zero bound states (ZBS). The basic process is shown in 
Fig. 3, where a d-wave S faces an insulator or a nor- 
mal metal N' along the (110) orientation (the a-b planes 
can be in principle rotated of an angle with respect the 
S/N interface different from 45°). An electron, coming 
for instance along the direction of the positive lobe of 
the order parameter, first suffers an ordinary scattering 
at the interface with a normal metal (N') and then is 
Andreev reflected towards the negative lobe of the order 
parameter. Then the hole will experience an ordinary 
scattering at the S/N' interface and will be reflected to- 
wards the positive lobe of the order parameter, where it 
will be Andreev reflected again. This process produces a 
constructive interference and as a consequence a resonant 
state, formally described by a pole in the Andreev ampli- 
tude Eq. ([!]). This manifests itself as a zero bias peak in. 
the density of states in S (zero bias anomaly) (ZBA) BED 

These arguments are essentially based on the retrore- 
flection property of Andreev reflection. On the other 
hand the constructive interference breaks down for high 
values of A/Ep: the electronic state created after two 
Andreev and two normal reflections will not propagate 
along the same trajectory as the initial one. As a conse- 
quence ZBS will be damped and the low voltage conduc- 
tance will be decreased. 

We have calculated the conductance of a tunnel NIS 
junction [Z\ — Z2 1) selfconsistently by the method 



decrease of the critical current density as a function of 
the tilt angle in HTS Josephson junctions, as experimen- 
tally observed. A damping of the resonances originating 
the zero bound states has been also found. 
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FIG. 4. The conductance for fixed value of Z\ — 0.5 and 
for different values of Z2 at temperature T = K. 

of RefM taking into account the angle mismatch of 
electrons and holes and the selfconsistent reduction of 
the pair potential at the interface. The angle averaging 
was performed by choosing a (5-function potential barrier 
with an angular dependence of the transmission coeffi- 
cient D{9) = D(0) cos 2 9. The results of calculations are 
presented in Fig. 3 for misorientation angle a = 45° and 
different values of the A/Ef ratio, giving some evidence 
of broadening of the ZBA. We also point out that in pj>p=i 
trast with ZBS broadening due to surface roughnesst§t2l 
the mechanism we propose has an intrinsic nature mainly 
controlled by the A/ Ep ratio. 

The problem of the ZBA has also been investigated 
in the case that Z\ is different from Z2. This difference 
is particularly meaningful in the formation of the ZBS, 
where both electron and hole scattering processes across 
the same interface are involved. In Fig. 4 we report the 
conductance in the regime A/Ep = for fixed value of 
Z\ = 0.5 and for different values of Zi at temperature 
T = K. We notice the appearance of peaks at finite 
voltages and the removal of the ZBA for some values of 
Z-2 > Z\. This means that the crossover from ZBA to 
bound states at finite voltages can also take place due 
to a different transparency of electrons and holes at the 
S/N interface. Such a crossover has been also predicted 
for ferromagnet -insulator - superconductor junctions by 
Kashiwaya et a/£j Our result is obviously independent 
of any ferromagnetic electrode or barrier and only relies 
on a possible barrier asymmetry for holes and electrons. 
Barrier asymmetry acts as a kind of filter which creates 
an electron-hole imbalance, thus reducing the probability 
of the formation of the ZBA. 

In conclusion, Andreev reflection has been theoreti- 
cally investigated in layered anisotropic normal metal / 
superconductor (N/S) systems in the case of an energy 
gap (A) in S not negligible with respect to the Fermi en- 
ergy (Ep). We have demonstrated that the combination 
of large gap and strong anisotropy leads to an intrinsic 
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